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Abstract 

Thermal inflation is an attractive solution to the cosmological moduli problem. However, domain 
walls may be formed after thermal inflation and some mechanisms are needed to eliminate the 
domain wall before it dominates the Universe. We point out that gravitational waves produced by 
the dynamics of domain walls may be observed by the pulsar timing experiments and future space- 
borne gravitational wave detectors, which provides a probe into the period of thermal inflation. We 
also show that the QCD instanton effect can effectively eliminate the domain walls with producing 
observable amount of gravitational waves. 



1 



I. INTRODUCTION 



While the supergravity and string theory are well-motivated candidates of the funda- 
mental theory, there often appear light scalar fields, called moduli, having Planck-scale 
suppressed interactions with matter. The moduli cause serious cosmological problems, since 
the lifetime is typically longer than 1 sec and the energy density stored in the form of coher- 
ent oscillation is too much. The decay of moduli usually injects significant amount of energy 
after the big-bang nucleosynthesis (BBN) begins, and it modifies subsequent cosmological 
scenarios. This is called the cosmological moduli problem Even if the moduli are 

heavy enough to decay before BBN, gravitinos produced by the moduli decay again may be 
problematic 

Thermal inflation {4- 7] is a short period of inflation which takes place well after the 
primordial inflation. It is regarded as a solution to the cosmological moduli problem, since 
the coherent oscillation of the moduli is sufficiently diluted by the exponential expansion 
during thermal inflation and subsequent entropy-production by the flaton decay. Thermal 
inflation is driven by a scalar field, called flaton, having a flat potential and trapped at 
the origin of the scalar potential due to thermal effects. This kind of scalar field can be 
embedded into supersymmetric (SUSY) theories. 

However, the most thermal inflation models suffer from a problematic domain wall (DW) 
formation after the phase transition in association with the end of thermal inflation 
This is because the flaton usually has Z n symmetry which is spontaneously broken at the 
true vacuum. Once formed, DWs eventually dominate the energy density of the Universe 
resulting in a cosmological disaster. Thus the Z n symmetry, which guarantees the flatness 
of the flaton potential, must be explicitly broken to some extent in order for the DWs to 
disappear before they come to dominate the Universe. Actually, an introduction of a small 
explicit breaking in the scalar potential can solve the DW problem jg, 9]. Such DWs can be 
a source of gravitational waves (GWs). 

In this paper, we study the GWs from the dynamics of DWs produced after thermal 
inflation. We show that the amount of the GWs may be within the reach of the pulsar 
timing experiments and future space-borne GW detectors. Thus the study of GWs provides 
a way to probe the period of thermal inflation. It is also shown that the Z n can be naturally 
broken by the QCD instanton effect in a simple class of thermal inflation model; in such a 
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model, the abundance of the GWs will be well within the observable range. 



II. THERMAL INFLATION MODEL 

Let us start with the model of thermal inflation. We introduce a gauge singlet superfield 
which takes the role of the flaton. 1 Imposing Z n symmetry, under which transforms as 
—7- e 2nt / n (j), we adopt the following superpotential 2 

W = + kcf>QQ + Wo. (1) 

Here, Q and Q are additional "quarks" which have quantum numbers for SU(3)qcd; they 
are thermalized when = 0, and give a thermal mass to 0. In addition, Wq is a constant 
which is related to the gravitino mass m 3 / 2 as Wo = m^^Mp (with Mp being the reduced 
Planck scale). In the following, in order to guarantee the flatness of the potential of 0, we 
take n > 4. We also comment here that, if Wo = 0, the superpotential has an i?-symmetry 
under which has a charge +2/n. Although such an i?-symmetry is explicitly broken by 
the constant term Wo, we will be interested in the case that the effect of the breaking is 
relatively small. Then, in the vacuum where ((f)) ^ 0, quasi Nambu-Goldstone boson, which 
we call i?-axion, shows up. 

Including the SUSY breaking effect, the zero temperature potential of the flaton takes 
the following form : 

V = Vo - + (» " 3) + ^) + ^ < 2 > 

If there are no sources for the A-parameter other than Wo, we have A = m 3 / 2 . The SUSY 
breaking mass parameter for is represented by —m 2 and is assumed to be negative. Here- 
after we take A real and positive. In addition, we take A m for simplicity. The vacuum 
expectation value (VEV) of is given by 

/mM-^V""" 21 



1 We use a same symbol for a superfield as its scalar component. 

2 Terms like <fi 2n , 3 ™ , . . . are also allowed, but they have little effects on the flaton dynamics discussed 
below and hence are neglected. 
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The flaton mass around the minimum is given by = y/2(n — 2)m, and Vq is determined 

as 

V = —m\\ (4) 
n — 1 

so that the cosmological constant vanishes at the vacuum. In the gravity-mediated SUSY 
breaking model, m is expected to be of the order of TeV. In the gauge-mediated SUSY 
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breaking model [10] , on the other hand, the scalar potential is slightly more complicated 
121 ] and the parameter dependence changes. Hereafter, we treat v and as free parameters 
without specifying the mechanism of SUSY breaking. 

Suppose that is trapped at the origin due to the large Hubble mass term during/after 
inflation. In such a period, Q and Q are massless and hence they participate to thermal 
bath. Thus the flaton obtains a thermal mass of ~ k 2 T 2 . As the cosmic temperature 
decreases, the energy density of the Universe is dominated by the potential energy of (p and 
thermal inflation begins. Thermal inflation ends at the temperature T en d ~ m for k ~ 0(1). 
Then the flaton oscillates around the minimum, and finally it decays into radiation. We 
parametrize the flaton decay rate as 

4:71 V 2 ^ 

with a numerical constant c, which depends on the decay mode. If decays into the Higgsino 
pair via the superpotential of k'<^H u H^ (where H u and are the up- and down-type Higgses, 
respectively) with k' ~ m^/v, then c ~ 0(1). On the contrary, c ~ O(10~ 4 ) for the decay 
into gluons via the loop effect. The reheating temperature after the flaton decay is estimated 



as 3 



90 n 1/4 



T R = y/TMp ~ 6.2 x 10 2 GeVVc 



228.75 



1/4 ( mj, y/2 pQiQGeV \ (g) 



VlTeV/ 



where g*(T) is the effective number of massless degrees of freedom at the temperature T. 
Notice that Tr is required to be higher than a few MeV for successful BBN. 

Let us estimate the modulus abundance in the presence of thermal inflation. We denote 
the modulus field by x an d its mass and energy density by m x and p x . Thermal inflation 



3 The flaton can also decay into the i?-axion pair. If the decay into the i?-axion is dominant, the final 
reheating completes by the decay of i?-axion. This situation is effectively incorporated in the parametriza- 
tion ([5]) by choosing the value of c appropriately In addition, the flatino, the fermionic component of the 
flaton, has a mass of y/n — lm, and hence the flaton cannot decay into the flatino pair. 
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efficiently dilutes the moduli through the late-time entropy production. The dilution factor 
is defined as the ratio of the entropy with and without the additional entropy production at 
the flaton decay It is given by (3, S] • 



30 V i* 

~'1.3xl0 15 



7T 2 #*( T end):T e 3 nd T R 



228.75 



g*(T ( 



end J 



/ !TeV \ / !GeV \ / v y 
V m J \ T R J Vl0 10 GeW ' 



The modulus abundance after the entropy production is 

pWprim) _ 1 1 / XQ N 

~) ~A8 R \W P 



s 



1X10 GeV ^(lTw) {-^—) [WGev {M- p ) 



(8) 



where Xo is the initial amplitude of the modulus. Here we have assumed that the modulus 
begins to oscillate before the inflaton decays and also that the inflaton decays before thermal 
inflation starts. Otherwise, the expression becomes more complicated 8|. We call this 
"primary" moduli for the reason discussed below. 

There is another contribution to the modulus oscillation [z]. Since the potential minimum 
of the moduli during thermal inflation may be displaced from the true minimum due to 
the Hubble mass correction, the secondary oscillation is induced after thermal inflation 
ends. The amplitude is estimated as 5\ ~ VoX*/ \ m \Mp). Hence the moduli abundance is 
estimated to be 

pN(sec) = m 2 x (5 X ) 2 /2 3T R 

S> " Vo 4 (9) 
-x 10 -^(^) 7/ X^) 2 («)fe) 2 - 

Here x* is the true minimum of the modulus while the origin is chosen so that the Hubble 
mass correction is given by H 2 \ 2 - This is called the "secondary" moduli. The final modulus 
abundance is the sum of the primary and secondary moduli. Which one dominates the 
moduli density depends on various parameters. BBN, cosmic microwave background, and 
diffuse X- and 7-ray background give stringent upper bound on the modulus abundance, 



depending on its mass. Comp 
effects can be found in Ref. 8]. 



ete analyses on the modulus abundance and its cosmological 
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III. DOMAIN WALL PROBLEM AND A SOLUTION 



In this section we discuss the problematic DW formation after thermal inflation, and a 
solution to the DW problem. 

After the decay, the flaton randomly falls into one of the n-degenerate minima. Since the 
VEV of spontaneously breaks the Z n -symmetry, DWs are formed after thermal inflation. 
The DW tension is estimated as 

a~— , (10) 

where 

^ = 2 (-- 3 )^^> (11) 



and 



■in 



2 



n{n — 3) 



Am (12) 



a 

gives the i?-axion mass. DWs obey the scaling low in which about one DW exists per Hubble 



horizon as far as the viscosity on the DW is negligible 13[. 4 Then DWs begin to dominate 
the Universe at 



#d om * ^j- (13) 



This is problematic unless a < (lMeV) 3 , since otherwise DWs dominate the Universe before 
the present epoch fl^ |. Thus there should be an explicit Z n breaking term which makes 
DWs unstable. 

Let us denote the Z n -breaking scalar potential by V e , whose possible origins will be 
discussed later. It generates a bias for the n minima of the original potential and the 
degeneracy among those minima are lifted completely. DWs are not absolutely stable under 
the existence of the bias, and eventually collapse when the bias energy density becomes 



comparable to the DW energy density 



151 ] . This happens when 



dec ^ (14) 

(7 



Therefore, in order for the DWs to collapse before they come to dominate (H^cc 3> -f^dom), 
we need the following condition, 



V. » ^. (15) 



4 Actually interactions of the DW with Q and Q cause friction, but this effect is found to be insignificant 
for the following analysis. 
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The DW problem is solved if this condition is satisfied. Now we describe possible origins of 
the bias term. 



A. Small explicit Z n breaking term 

Let us introduce the following superpotential 

in addition to flTJ), where n and i are relatively prime numbers and e is a small coefficient. 
In order for this term not to change the flaton potential significantly, we require 

v e v n I ( v \ n ~ e j . 

« 77777^ e<- T7 • ( 17 ) 



£M l ~ 3 nM n ~ 3 n\M 
Then the following additional scalar potentials are generated, 

±t ±n-\X*l-\ 

v * = ^ ~ 3 )^=3 + + h - c - ( 18 ) 

The second term dominates for m > A and this yields the bias for the n minima of the 
original potential as 

n-t-2 

The original Z n symmetry is explicitly violated by the bias, and hence the degeneracy of the 
n minima of the scalar potential is broken. The condition (Tl5|) is rewritten as 

A v 2 ( v \ n ~ e 

e>> ^w P U ■ < 20 > 

If this condition is satisfied, DWs collapse before they come to dominate the Universe and 
there is no DW problem. 



B. QCD instanton effect 

Rather simple scenario to break the Z n symmetry is to use the QCD instanton effects [Tj 
This is economical in the sense that the term <pQQ in (JTJ, which is needed for giving rise 
to a thermal mass for the flaton field, also works as a source of the bias for eliminating the 
DWs. 5 



This solution may be incompatible with the Peccei-Quinn mechanism [3| for solving the strong CP 



problem Yjj, |18(. In that case, we may arrange the model so that the heavy quarks Q(Q) have a charge 



of hidden QCD and its scale is around GcV-TcV. 
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The Z n symmetry in ([!]) has an anomaly for SU(3)qcd- Consequently, the QCD instanton 
effect works as a bias, which lifts the classical degeneracy of the n vacua. The bias due to 
this effect is estimated to be [l^j 

(m u + m d y 

for temperature below Aqcd ~ 200MeV, where f n and are the pion decay constant and 
mass, and m u (md) is the current mass of the up (down) quark. At higher temperature, this 
effect is suppressed by high powers of rfl^ence the h iaS due to the QCD Canton 
effect is only turned on at T ~ Aq CD 
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We need two conditions for this mechanism to work successfully. One is Eq. ( fl5l) for 
DWs collapse due to the bias. The other is H& cc < Aq CD /(T^Mp), since otherwise the QCD 
instanton effect is highly suppressed and it is not suitable for the bias. 6 These two conditions 
are written as 

TlM P <a< A 2 QCD M P (22) 

for Tr < Aqcd- Therefore, for example, M ~ 10 15 GeV and A = 7713/2 ~ m ~ lGeV are 
good choices for n = 4. In this choice we have H^ ec / H^ om ~ 0(1). As will be seen, this 
predicts observable GWs. 

It should be noticed that M cannot take an arbitrary value in order to solve the cos- 
mological moduli problem. As seen in Eqs. (IE|) and ([9]), the moduli abundance crucially 
depends on v, which is determined by the cutoff scale M. Comparing with the BBN bound 
Px/ S ~ 10~ 14 GeV for the moduli of lOGeV < m x < ITeV 
may be allowed for moderate parameter choices. 



M > 10 15 GeV for n 



IV. GRAVITATIONAL WAVES FROM COLLAPSING DOMAIN WALLS 

Now we are at the position to discuss the main subject of this paper, which is the GW 
production from DWs. We will show that the GWs from the DWs after thermal inflation 
may be observable by future experiments. 

The DW network has complicated structure. Large amount of energy is localized around 
DWs, which becomes the source GWs. As a result, considerable amount of GWs can be 

6 The latter condition is based on the assumption that the Universe is flaton dominated at the QCD phase 
transition. If the Universe is radiation dominant at the QCD phase transition, the condition is severer. 
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emitted by the relativistic motion of DWs. Frequencies of the GWs correspond to typical 
scales of the DW motion, which ranges from the Hubble scale to the size of DW width. 
Detailed calculations of the GW spectrum from DWs using three dimensional lattice 



simulation were performed in Refs. 24H26j . According to recent studies 25|, |26] , GWs have 
rather broad spectrum which extends from the horizon size at the DW collapse for lower 
frequency side to the DW width for the higher frequency side. The GW spectrum, in terms 
of Ogw(/) = (dpcw / dlia f) / Pco (where pcw(f) is the GW energy density with frequency 
/ measured at present and p c o is the present critical energy density) is almost flat between 
these frequencies. The energy density of gravitational waves at the collapse of DWs can be 
estimated as 



Pcw(H dec ) ~ G 



/v 



-^dec -^dec 



a 



(23) 



where Gn is the Newton constant. Here we have substituted Mow = a H A ^ c . Then, Qqw is 



given by 



^GW(/) - 
9* (f dec 



dec J 



9*0 



9*0 
9*s0 



9*s0 



_g*s(T ( 



dec ) 



4/3 



H, 



dom 



H, 



dec 



4/3 



H, 



dom 



R 



2/3 



_9*s (-^dec 

for / ed g C < f < /peak defined below, where O 
with three massless neutrino species density. Numerically we obtain 

2 



for H R >H dcc , (24) 
for H R <H dcc , (25) 



H dcc J \ H dcc 

= 8.5 x 10~ 5 is the present radiation energy 



n GW (f) ~ 2 x io- 
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" 228.75 " 


1/3 


H dom \ 




.^(^dec). 


1 


V Hdcc J 



f2 GW (/)-2xRr 



228.75 



9*(T C 



dec ) 



1/3 



H, 



dom 



H, 



dec 



2/3 



for H R >H dec , (26) 
for H R <H dec , (27) 



for /edge < f < /peak- The typical frequency at the lower side, corresponding to the horizon 
size at the DW collapse, called the "edge" frequency (/edge) in Ref. [25] . is given by 

-f^dec a(H dec ) 



/edge 

which is estimated as 

/edge ^ 3 X 10~ 5 HZ 

" (rp NT 1/6 

9*{Tr^ 



2tt a Q 



(28) 



g*(T dec ) 


1/6 

I 


< Tdec \ 


_ 228.75 _ 




v lTeVy 



/edge ^ 3 X 10" 5 HZ 



228.75 



7i 



R 



ITeV 
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H, 



dec 



H 



1/3 



for H R > H dec , 
for H R < if dec, 



(29) 
(30) 



where T dcc is the temperature at which the DWs collapse. On the other hand, the "peak" 
frequency (/peak) corresponding to the DW width at the collapse redshifted to the 
present time, is given by 



peak 



w 1 a(H dec ) 



where w 



2n a 

~ m~ x is the DW width. This is estimated as 



(31) 



a 



/ peak ~ 1 x 10 iU Hz 



228.75 



/ peak ~ 1 x 10 iU Hz 



228.75 



_g*(Tdec) 

( m„ 



1/3 



(■ 



1/3 



VlTeV 



VlTeV 

m 



/lTeV\ 



T, 



dec 



H 



R 



H, 



dec 



2/3 



for H R >H dcc , (32) 
for H R <H dcc . (33) 



We approximate that the GW spectrum is flat between / e d ge and / pC ak, although the simu- 
lation shows a bit preference for the blue spectrum 25j. In that sense, our estimate of the 
GW is conservative. It may be within the observable range of GW detectors such as pulsar 
timing array or space laser interferometers 



Fig. [Tj shows the edge frequency / c d gc (left) and Oqw (right) as a function of H dec / H dom , 
which parameterizes the epoch of the collapse of DWs. Here we have taken m = ITeV, 
A = lOOGeV, c = 1, M = Mp and n = 4,5,6. It is clear that the GW energy density 
becomes larger for DWs that collapse closer to the epoch of DW domination. 

Fig. [2] shows contours of Og W , / e d ge and T R for n = 4 (top) and n = 5 (bottom). In 
the gray region denoted by "no DW", the bias is so large, i.e. ifdec > v^Vo /Mp, that no 
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3 11 ] and 



DW formation is expected. The current limits from pulsar timing experiments 
the LIGO experiment do not give stringent constraints on the model. In the same 
figure, we also show the prospects of the discovery reaches of up-comming projects. Projects 
such as SKA 



331 ] wil 



the advanced LIGO 



_probe some parameter regions. Ground based GW detectors such as 



34j and LCGT 35j are also sensitive to the higher frequency GWs of 



/ > 100 Hz. 
such as LISA 



3ere a correlation analysis of 1 year is assumed. Space based GW detector 
will also be sensitive to the GWs from DW collapse for wider parameter 



37j. DECIGO |38] may cover the whole parameter 



regions. Sensitivities are found in Ref. 
region with Qqw ^ 10~ 18 . The scenario with the QCD instanton as a bias term corresponds 



to the lower-left edge of Fig. |2] and hence future pulsar timing and/or space-borne GW 
experiments will be able to detect GWs. 
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'dec' dom n dec' dom 



FIG. 1: The typical GW frequency /edge (left) and Oqw (right) as a function of H& ec / -ffdom- We 
have taken m = ITeV, A = lOOGeV, c = 1, M = M P and n = 4, 5, 6. 



V. DISCUSSION AND CONCLUSIONS 



Some comments are in order. It was proposed that the first order phase transition trig 
the end of thermal inflation, and correspondingly bubbles of true vacuum are created 



gers 



39] 



Bubbles expand into the surroundin g fa lse vacuum regions and collide with each other, 



generating GWs with observable level [39j . This predicts higher frequency GWs than studied 
in this paper and hence will provide further information on the thermal inflation model. 

We also comment on possible origin of the baryon asymmetry in the model. Since thermal 
inflation dilutes away the pre-existing baryon asymmetry, as well as the moduli, we need 

lis regard, a 



41 



45]. 



some mechanisms to create baryon asymmetry after thermal inflation. In t 
variant type of the Affleck-Dine mechanism 40[ works after thermal inflation 

One might notice that the heavy quraks Q and Q may be stable and cosmologically rele- 
vant 46|]. Since they are once in thermal equilibrium during thermal inflation, its abundance 



just after the flaton decay is given by 



Yn 



s 



end R 

m 2 v 2 



10 



-17 



R 



( m <t> \ 
UTeV/ VlGeV 



10 10 GeVV 



(34) 



Even if it is less than the DM abundance, BBN constraints on the strongly-interacting relic 



particles are much more stringent 



481] . 7 This constraint can easily be evaded if Q and Q 



' The enhancement of the annihilation cross section due to the i?-hadron formation (^J may help the 
situation if the flaton decays after the QCD phase transition. 
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M [GeV] 

^22 




r R = lOOGeV 



T R = ITeV 



'"'dec^dom 



M [GeV] 

.22 




T R = lOOMeV 



T R = lOGeV 



'"'dec^dom 



FIG. 2: Contours of Ogw (solid line), /edge (dotted line) and Tr (dashed line) for n = 4 (top) and 
n = 5 (bottom). The regions explored by LISA, SKA and advanced LIGO are also shown. DECIGO 
may cover the whole parameter region with ficw ^ 10 -18 . Here we have taken m = ITeV, 
A = lOOGeV and c = 1. 
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are in complete multiplets of the GUT group, which is natural from the viewpoint of gauge 
coupling unification. If Q and Q are embedded into fundamental and anti-fundamental repre- 
sentation of SU(5)gut, for example, they should be acompanied by non-colored components 
(denoted by L and L). The renormalization group effect makes the colored components 
heavier than non-colored ones at the low energy. Then, Q and Q decay into L and L via 
dimension five operators suppressed by the GUT scale. The decay rate is estimated to be 
T ~ m g/^GUTi where M GUT ~ 10 16 GeV is the GUT scale. The decay rate is large enough so 
that colored components decay before BBN. Similarly, L splits into the electrically charged 
and neutral components, where the former is heavier and decays into the latter. Notice that 
the neutral component may be a candidate of superheavy DM. 8 

Let us consider how robust the DW formation after thermal inflation is. In the simplest 
model, the flatness of the flaton potential is ensured by a Z n symmetry as in ([I]). In this 
case DWs necessarily appear. On the other hand, one can impose a global U(l) symmetry 
to make the flaton potential flat. The flaton potential is stabilized by introducing another 



singlet scalar wit 
flaton soft mass 



i some U 



45 



51 



1) charge or by taking account of the radiative correction to the 
55| . Then, the spontaneous breaking of the U(l) symmetry due to 
the VEV of the flaton results in the formation of cosmic strings instead of DWs. Even so, the 
global U(l) symmetry may be anomalous and may not be a good symmetry at the quantum 
level. If only a discrete subgroup Zn of the U(l) symmetry remains at the quantum level, 



DWs still exist 



. In particular, if N > 2, the DW problem still exists. On the contrary, 



if N = 1, DWs are bounded by strings and they disappear just after the formation 57H61|. 



It is also possible that the U(l), which guarantees the flatness of the potential, is a gauge 
symmetry 6j. The flaton takes the role of the Higgs field that gives a mass for the U(l) 
gauge boson. Then, the U(l) should be an exact symmetry and DW does not exist. Instead, 
cosmic strings are formed which are less harmful than DWs. Cosmic strings formed after 



thermal inflation has a width of the TeV scale, while the tension can be much larger 



kind of thick strings have characteristic cosmological implications as studied in Refs. 62 



'his 



63] 



Finally, we comment on another possible source of GWs, which is the DW related to the 
breaking of Peccei-Quinn symmetry [l^] . In order to avoid the DW problem, it is often the 



The constraints from direct detection experiments can be avoided if the mass of the neutral component 
is heavier than ~ 10 5 GeV [50 1. 
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case that axion models with N = 1 is considered. In such a case, DWs are surrounded by 
cosmic strings as we have mentioned; then, DWs disappear just after the production. This 
happens when T ~ Aqcd, so we expect that the GWs with the frequency of / ~ 10~ 8 Hz are 
emitted from the dynamics of axionic DWs. We have estimated the amount of such GWs, 
and found that f^ow is a few order of magnitude smaller than the sensitivity of the future 
pulsar timing experiments. 

To summarize, a class of thermal inflation models is inevitably associated with DW 
formation, which reintroduces a cosmological disaster. Thus we need a mechanism to make 
DWs unstable. If the DW energy density at the epoch of its collapse is large enough, GWs 
produced by the DW collapse may be observable in future experiments. We have pointed out 
that such a mechanism is naturally built-in in the simplest model where the flaton couples 
to extra quarks as in (TjQ); in such a model, the Z n symmetry has an anomaly for SU(3)qcd 
and the QCD instanton effects serve as a bias for lifting the degeneracy among n minima. In 
this case the GW amplitude becomes so large that it can be well within the reach of future 
experiments. 
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